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| Abstract. In the presence of the primordial magnetic held, initial vector (vorticity) 

perturbations produce cosmological Alfven waves and leave imprints on cosmic 
microwave background (CMB) temperature and polarization anisotropy. 

We have investigated imprints of cosmological Alfven waves in CMB anisotropy. For 
data constraints, we have used the power spectrum of the recent CMB observations, 
and correlations estimated from WMAP Internal Linear Combination (ILC) maps. 
Our analysis shows 3cr evidence of cosmological Alfven waves. Using the 3cr limit from 
our analysis and the Alfven velocity limit from the total energy density constraint, we 
impose a lower bound on the amplitude of primordial vector perturbation: 4 x 10 -12 
at fc = 0.002/Mpc. 
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O 1- Introduction 

G\ ' 

There are strong observational evidences on the existence of large-scale magnetic fields, 
whose strength is on the order of micro Gauss [H El El Hj. These magnetic fields are 
believed to be seeded by a small Primordial Magnetic Field (PMF) [SI El 13 IE]- Primordial 
Magnetic Field may be considered to consist of a homogeneous and inhomogeneous 
(stochastic) magnetic fields. In this paper, we focus on the effect of Primordial 
Magnetic Field on CMB anisotropy associated with vector perturbations. Primordial 
vector perturbation, whose source includes topological defects and inhomogeneous 
primordial magnetic fields 0, EUS [11], decays rapidly with the expansion of the 
Universe, and therefore, does not leave observable signature on CMB anisotropy [T2"l[T3] . 
However, primordial vector perturbation in the presence of a homogeneous primordial 
magnetic field induces Alfven waves in the cosmic plasma, which may leave observable 
imprints on CMB anisotropy via a Doppler effect and an integrated Sachs- Wolfe effect 
[EU EES! ESI EE3 EES EES]- To be specific, cosmological Alfven waves create correlations 
between cti m , ai±Ai, m ±Am, where Al = 0, 1, 2 and Am = 0, 1, 2 [20]. Over the past several 
years, there have been great success in measurement of Cosmic Microwave Background 
(CMB) anisotropy by ground and satellite observations [21l [22l [23l [2^1 [251 [26l [2TI [28], 
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and observational imprints of cosmological Alfven waves in CMB anisotropy have been 
studied by several authors [HI [201 [29], [301 EE]- By investigating the effect on power 
spectrum and signature correlations in the recent CMB data, we attempt to constrain 
cosmological Alfven waves, and also impose a lower bound on the primordial vector 
perturbation. 

The outline of this paper is as follows. In Section [21 we briefly review the effect 
of cosmological Alfven waves on CMB anisotropy. In Section [3J we discuss statistical 
properties of CMB in the presence of the Alfven waves. In Section HI we build the 
statistics sensitive to the imprints of Alfven waves. In Section [51 we analyze the recent 
CMB data and present the result. In Section El we make summary and conclusion. In 



Appendix A, we discuss the effect of incomplete sky coverage on the analysis of Alfven 



wave imprints. 

2. The effect of Alfven waves on CMB anisotropy 



According to linearized Einstein equations [32] , vector perturbation decays rapidly with 
the expansion of the Universe, and therefore, does not leave observable signature on 
CMB anisotropy, unless the initial values of vector perturbation are unusually high 
[T2"l [T3] . In the presence of a homogeneous magnetic field, vector perturbation induces 
Alfven waves in cosmic plasma, and leaves observable imprints on CMB anisotropy 
via a Doppler effect and an integrated Sachs- Wolfe effect [18J. Durrer, Kahniashvili 
and Yates [H] (hereafter, DKY) showed that cosmological Alfven waves generate the 
fractional CMB anisotropy for a Fourier mode k: 
ST 

— (n,k) « n • ft(k,r] dec ) = n • O v A kr] dec B • k, (1) 

where n and B denote sky direction and a homogeneous magnetic field direction 
respectively, f2(k, T]d ec ) is the Gauge invariant linear combination associated with 
vector perturbations, rjdec denotes the conformal time at the moment of baryon-photon 
decoupling, va is Alfven wave velocity and To is the CMB monopole temperature 
2.725[K] [33]. DKY assumed that vector perturbations are initially created by some 
random process and have the following statistical properties over an ensemble of 
universes: 

(O J (k)Q J (k)) = (^-^%)P(A;), (2) 

P(k), which is the power spectrum, is assumed to follow a simple power law: 

k n 

p (k) = A v —^, (3) 

Av 



where ko is a pivot wavenumber and we set it to 0.002/Mpc, which is equal to the WMAP 
team's pivot wavenumber for scalar perturbation [31]. Possible sources of primordial 
vector perturbation include inhomogeneous stochastic primordial magnetic fields and 
topological defects 0CEO1QI]. 
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3. CMB statistical properties 

The temperature anisotropy 5T(6, <p) is conveniently decomposed in terms of spherical 
harmonics Y im (9, (j>) : 



5T(9,<fr) = ^a /m F Zm ( 



(4) 



lin 



where a/ m are the coefficients of decomposition, and 9 and <fi are a polar and an azimuthal 
angle. Kahniashvili, Lavrelashvili and Ratra [20] (hereafter KLR) showed that the CMB 
anisotropy in the presence of cosmological Alfven waves has the following statistical 
properties: 

1(1 + 1) 



(a* lm a lm ) = + (3cos 2 9 B - I) 



(21 - l)(2/ + 3) 



x 



1(1 + 1) + (I 2 + I - 3) cos 2 9 B 2 

m 



( a *m a l,m+l) : 
( a *lm a l,m+2) '■ 
{0-Lm a l+2,m) 



3 COS 2 9 B — 1 



sin 29 B exp[—i<j) B ]I 1 / d™' m+1 , 



\ sin 2 9 B exp[-2i(f) B )I 1 / d^ ,m+2 , 



1(1 + 1) 



rU 



(3 cos 2 9 B - l)I l / +2 <% 



l+2i 



(a ltm ai + 2, m +i) = sm29 B exp[-i(pB\Ia d lt ' +2 , 
(al m ai +2 , m -i) = sin29 B exp{i(f) B ]I l / +2 d™'™' 1 , 

(a lym a l+%m+2 ) = - -sin 9 B exp[-2i<j) B }I^ d lt > +2 , 



(0'l m ai + 2,m~2) 



• 2 n [n-i i tUI+2 im,m—2 

sin 9 B exp[2z(pB\I c l a n + 2 , 



(5) 

(6) 

(7) 

(8) 
(9) 
(10) 

(11) 
(12) 



where Cf is the power spectrum in the absence of Alfven waves, 9b and <pB are the 
spherical coordinate of a PMF direction B, and o?^' m (|m| < I, \m'\ < V) are 

— - — ^rr-r— — r (m -\ — ) \/ (I — m)(l + m + 1), 
(21- 1)(2/ + 3) V 2 ;v v A ; ' 



jm,m+l 
a Ll 



a i,l+2 



im,m+l 



tm,m—\ 
a l,l+2 



P + 1-3 

(2/ — 1)(2Z + 3) 
(Z + 3)Z 



2(2/ + 3)v/(2/ + l)(2/ + 5) 



^/(Z - - m - 1)(Z + m + 1)(/ + m + 2), 

y/((l + l) 2 - m?)(l - m + 2)(l + m + 2), 



(1 + 3)1 



2(21 + 3)^(21 + l)(2/ + 5) 



(Z + 3)Z 



2(2/ + 3)v/(2/ + l)(2/ + 5) 



VT((/ + l) 2 - m 2 )(l + m + 2)(l + m + 3), 
V((Z + l) 2 - m 2 )(Z - m + 2)(l - m + 3), 
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+2 (1 + 3)1 /(Z + m + 4)! 

a Ll+2 



2(2/ + 3)^(2/+ l)(2/ + 5) V (Z + m) 



and 



d m,m-2 = (* + 3)Z / (Z -771 + 4)1 

M+2 " 2(2/ + 3)v/(2/ + l)(2Z + 5) V (Z-m)! 



7 «< = HZo y fc2 (P Q ^ exp (-2-|) (^)^MiKM*, 

x nv+2 exp ( - lu ji(x)j[{x) dx. 



2T 2 P ^ f Vdec y /".+»/ 2x 2 



7T (/co^o)"^ 3 \ Vo J J \ (k D Vo) 

(13) 

?7o is the present conformal time, and is the Alfven velocity, which is given by 

v\ = - B » , VA ~4xlO- 4 -— § , (14) 

A A-K(p r +p r y 10" 9 Gauss v 1 

where B is the magnitude of a homogeneous PMF, and p r and p r are the density and 
the pressure of photons |18j . ho denotes the comoving wavenumber of the dissipation 
scale, due to photon viscosity and given by ~ 10 jr}dec [IB]- While DKY and KLR have 
obtained analytic results by neglecting the damping factor exp (— 2k 2 /k^), we have 
computed Eq. [13] numerically. The damping effect is getting significant on multipoles 
I > 500 [IB]. 

Primordial magnetic fields affect the expansion dynamics of the Universe, because 
of anisotropic pressures associated with primordial magnetic fields [35]. Considering 
the shear anisotropy due to primordial magnetic fields and the observed magnitude 
of CMB anisotropy (~ 10~ 5 ), Barrow et al. (1997) have derived an upper limit 
B < 2.27 x 10~ 9 h/75 Gauss on a homogeneous primordial magnetic field [35J. Using 
this upper limit and Eq. [TJJ DKY showed B < 10~ 8 Gauss and v A < 10~ 3 [15] . 

It is worth to note that correlations are invariant under the parity inversion of a 
PMF direction B. In other words, (6b,<Pb) — 6b,4>b + ^)) does not affect Eq. El [71 
El El dOl EH and [T2l note that it is the degeneracy of correlation, not that there should 
be two physical directions. 

4. Estimators 

In order to constrain cosmological Alfven waves, we construct the following statistics: 
Ct =(2l + l)- 1 Y J oT lm a lm , (15) 

m 

D l = Yl a *lm a l,™+l, (16) 
m>0 

Df = (2/-l)- 1 ^a i >,, m+2 , (17) 
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Df = (2/ + l)- 1 ^a* m a, +2 , m , (18) 

m 

Df = (2/)- 1 ^a* m a z+2 , m+1 , (19) 

m 

D? = (2/-l)- 1 ^< m a i+2 , m+2 . (20) 

m 

Because of the reality condition a;_ m = (— l) m ai m *, a* m a/ >m+ i is always zero. 

m=-2l 

Therefore, we defined Dj to sum only terms of m > 0. There also exist the following 
indentities because of the reality condition: 



in 



Q>l+2,m+l — — 



2^ a t' m a ' + 2 : m " 1 

m 



«Z+2,m+2 



E 



a Lm a l+2,m-2 



(21) 



(22) 



Therefore, Y^ m a tm a i+2,m-i and Y^ m a *i,m a i+2,m-2 are redundant and excluded from 
consideration. Note that the statistic used in [181 EH] is equivalent to Eq. [181 Using Eq. 
El El El M and [HI we may easily show that the expectation values of C\ and D\ are 



C, 



A 1 



A 2 



D 



c o + K[±i) /W 



d • 



sin2^ B exp[— i<p B ] -j- d 



>m,m+l 

m>0 



- sin 2 B exp[-2z</> Bj ^ 



T i,i 



•m,m+2 



',1+2 



(3 cos 2 6 B - 1) 



2/ + 

rU+2 



Z?f = sin2#e exp[— 



2/ 



A 5 



^sin 2 B exp[-2i0 Bj o 



Eim,m 
a i,l+2i 

m 

E< 

r, 
/ 



Z+2 J 



m 

U+2 



iE<* 



m,m+2 
1,1+2 ■ 



We may also show that the variance of D\ is 
Var(Ci) ~2/(2/ + l)(Cf + iV i ) 2 , 

r 1 (c? + jv / ) a , 



Var( J D/) 
Var(A) 

Var(L>f) » (2/ + l)- 1 ^ + iV0(C° +2 + iV i+2 ), 
Var(Df) 
Var(A) 



(23) 
(24) 

(25) 
(26) 
(27) 
(28) 

(29) 
(30) 
(31) 
(32) 
(33) 
(34) 



(2l)- 1 (Cf + N l )(Cl 2 + N l+2 ), 
(2/-l)- 1 (Q° + iV0(C° +2 + iV, +2 ), 

where Ni is noise power spectrum, and we have neglected correlations between distinct 
spherical harmonic modes (i.e. (a* m ai' m i) ~ Cf Su'S mm i). 
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Though the underlying distribution for primordial vector perturbations is not 
necessarily Gaussian, the distribution function of D\ — D\ tends to Gaussian by the 
central limit theorem [36]. Hence, the likelihood function of D\, given the data D\, is 

C{D\{X a )\D\) = 1 ex P [-i(D - D(A Q ))t M" 1 (D - D(A a ))]. 
(2tt)t|M|2 2 

(35) 

M is a 5 / max x 5 Z max nearly diagonal matrix, whose diagonal elements correspond to 
the variances of Dj, . . . , Df: M m diag (Var(D/ m . n ), . . . , Var(Df max )) . N is the data 
vector size (i.e. 5Z max ), and X a are the parameters associated with Alfven waves : 
A Q G {A 

v v \i n vi @Bi 0b}- F ot 1° w l-i the deviation of Dj — Dj from Gaussian distribution 
is getting non-negligible. However, it does not affect the best-fit values, but only the 
error bars, since M is nearly diagonal, and thereofore the best-fit parameters are, in fact, 
determined by minimizing \ 2 - Since ((C/ — C\) [D\ — D\)) is negligible in comparison to 
the variance of C\ or D\, we may construct a full likelihood function as follows: 

CiQlQ) x £(Dl(\ a )\Dl), (36) 

where C{Ci\Ci) is the likelihood function associated with CMB power spectra, whose 
approximate and optimal expression is found in [37]. Note that cosmological Alfven 
waves affect CMB power spectrum as well, and hence Eq. [23] should be used for 0% in 
estimating C{Ci\C{). Since CMB power spectra and correlations (i.e. C\ and D\) possess 
non-trivial dependence on cosmological parameters (i.e. standard parameters + Alfven 
wave parameters), we will resort to Markov Chains Montel Carlo (MCMC) likelihood 
analysis [374 EH] • The discussion on MCMC analysis and dataset used will be give in 
the next section. 



5. Analysis of the recent CMB observation data 

We have investigated the effect of the WMAP team's KQ75 and KQ85 mask on 
correlations (i.e. D\ ), and found that significant amount of fictitious correlations are 



produced by the masks (refer to Appendix A for details). Therefore, we find sky maps 



which require a Galactic mask, are not suitable for our analysis. Therefore, we have 
restricted our correlation estimation to whole sky maps obtained by variants of Internal 
Linear Combination (ILC) method: The WMAP team's ILC map (WILC), Harmonic 
ILC (HILC) and Needlet ILC (NILC) [391 US El S2J 113] . In Fig. HJ we show the 
temperature power spectra of the ILC maps with that of the WMAP concordance model 
[3H E3]. Noting that the power spectra of ILC maps are in good agreement with the 
model up to I < 200, we have computed D\ up to / < 200. We have also made the 
analysis with D\ <2b0 for HILC, and D\ <2b0 and -D;< 30 o respectively for NILC. We have 
obtained the consistent results with the D\ <2m analysis. The confidence intervals are 
smaller, but the best-fit values turn out to be similar. From now on, all the results are 
from the D\ <2m analysis. 
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i 

Figure 1. the power spectra of ILC maps and the WMAP concordance model 

Table 1. cosmological parameters of ACDM + (A v v\, n v , 9b, 4>b) constrained by 
power spectra (WMAP5YR + ACBAR + QUaD) + correlation (ILC) 



parameter 


WILC 


HILC 


NILC 


Q b h 2 


n no'i+0-002 

U.UZO_ .0Q2 


n n24 +0 ' 002 


0.024jS 


n c h 2 


0.109j£°£ 


n in4+o.o3 

U.1U^_q 01g 


O.lOSjffi 


T 


0.089+°°,f 




0.089±g-gS 


n s 


i-oi5±8:2g 


1.028i^| 


l.OlSj^ 1 


dn s /d\nk 


-0.0081H1 


-o.oo8ig:8g 


u.uuy_ 108 


log[10 10 A s ] 


3 . 03 X q 2y 


9 O7c;+0.228 
' d -0.234 


3.026jSig 


A sz 


1-631S 


1 87+0.133 
L -° ' -1.867 


1 r 7+O.43 


H [km/s/Mpc] 


74.37+^.65 




74.46l^;|I 


10 9 A v v\ 


q -i q+5.66 
0. 10_2 42 


4.161*3 


3.87±i"?f 


n v 


-4-35 + °f 


^■^ d -0.71 


-4.49+ ,;? 3 


e + B 


1 c 70 + 25 
i0 -' -15.7 


1 fi o+31.7 
L - u -1.6 


1 1 70+27.7 
lL -< -11.7 


4>i 


1 1 90+247.9 
iiz -111.8 


1 1 90+246.6 
HO. Z _xi3.l 


1 1 qo+240.9 
iiJ -119 


Ob 


164.3°i^- 7 


178.4°1^ 6 7 


168.3°i^ 


<Pb 


9Q90+247.9 
zyz _ 1118 


9QO 90+246.6 

zyo.z _ 1131 


2 qqo+240.9 
zyy -119 



By making small modifications to the CosmoMC package [38J, we have included the 
parameter X a G {A v v\,n v ,9B,4>B}, and computed the likelihood function given by Eq. 
1361 Note that the power spectrum as well as correlations depend on A v v\ and n v (see 
Eq. [23l) . For data constraint, we have used the CMB power spectra of the recent CMB 
observations (WMAP5YR + ACBAR + QUaD) [21 E21 ES EEJ EEJ E3 EE], and the 
correlations estimated from ILC maps. In Table [H we summarize la constraint and 
the best-fit values. As discussed in Sec. IU our estimators are insensitive to the parity 
of a PMF direction. Hence, we quote two best-fit PMF directions, which are equally 
likely. In Fig. El we show the marginalized likelihood (solid lines) and mean likelihood 
(dotted lines) of {A v v\, n v , 9b, <Pb}- In Fig. El we show likelihood in the plane of A v v\ 



Cosmological Alfven waves and the observational bound on the vector perturbation 8 





Figure 2. Likelihood of Alfven wave parameters: WILC, HILC and NILC (from the 
left to the right), normalized to its peak. 



versus n v , which are highly correlated Alfven wave parameters. We find there exists 
little correlation between Alfven wave parameters and ACDM parameters, except for 
the scalar perturbation amplitude A s . In Fig. [3J we show marginalized likelihood in the 
plane of A v v\ versus log[10 10 A s ]. Marginalized likelihood of PMF direction is shown 
in Fig. HI Note that a pair of directions associated by (9b,4>b) ^ (tt — 9 b ,4>b + 
possess equal likelihood. As shown in Table Q] and Fig. [2, [31 HI The result from various 
ILC maps are similar with each other. Since each ILC map contains distinct residual 
foregrounds and noise, consistency among the analysis with various ILC maps indicates 
that our results are affected insignificantly by residual foregrounds and noise in ILC 
maps. 

For good exploration on tails of the parameter A v v\, we have also set CosmoMC to 
check convergence for confidence intervals of A v v\. In Table |2l we show 3cr constraints 
on A v v\. Using our 3a limit on A v v\ and the upper limit on va (refer to Sec. 3), we 
find 5 x 10~ 4 < A v (refer to Eq. [3] for the definition of A v ). Taking into account Eq. 
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Figure 3. Marginalized likelihood in the plane of A v v\ vs {n v , log[10 10 A s ] }: WILC, 
HILC and NILC (from left to right), solid curves denote la and 2a contours. 




Figure 4. Marginalized likelihood of PMF direction: WILC, HILC, and NILC (from 
left to right), normalized to its peak 



|3l we find the amplitude of primordial vector perturbation is equal to A v kl. Hence, we 
impose the lower bound 4 x 1CT 12 on the amplitude of primordial vector perturbation. 



Table 2. A v v\ estimation at 'da confidence level 



Power Spectra + ILC 




3(7 


confidence interval 


WILC 


5.8 x 


io- 


10 < A v 


v\ < 1.01 x 10~ 8 


HILC 


5.73 x 


10' 


- 10 < A v 


v\ < 0.99 x 10~ 8 


NILC 


5.05 x 


io- 


- 10 < A v 


v\ < 1.09 x 10- 8 



If there was a primordial magnetic field at the recombination epoch or after the 
Universe was re- ionized, it would induce Faraday rotation in CMB polarization pHl |4"5] . 
Forecast on the PLANCK data constraint shows that we will be able to constrain a 
primordial magnetic field ~ 2 x 10~ 9 Gauss by investigating Faraday rotation effect in 
polarization data [4"5] . However, we have not considered Faraday rotation effect in this 
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work, since the signal-to-noise ratio in the currently available CMB polarization data is 
not high. 

6. Discussion 

We have constrained primordial magnetic field and primordial vector perturbation, by 
investigating the imprints of cosmological Alfven waves on CMB anisotropy. We find 
there is little degeneracy between cosmological Alfven wave parameters and ACDM 
parameters except for a scalar perturbation amplitude A s . The results obtained with 
various ILC maps are consistent with each other. Therefore, we believe our results have 
not been affected significantly by residual foregrounds or anisotropic noise in ILC maps. 
Using our result and the upper bound on va from the total energy density constraint, 
we impose a lower bound on the primordial vector perturbation 5 x 10~ 4 < A v at the 
pivot scale k = 0.002/Mpc. 

Since Alfven velocity v a is proportional to the magnitude of a primordial magnetic 
field, direct constraints on Alfven velocity va may be imposed by investigating Faraday 
rotation effect in CMB polarization data [44j [45] . When the data from the PLANCK 
surveyor [16] are available, we are going to constrain Alfven waves and Faraday 
rotation simultaneously, and hence impose a stronger constraint on the primordial vector 
perturbation. 
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Appendix A. incomplete sky coverage 

Incomplete sky coverage produces fictitious correlations by destroying the orthogonality 
of spherical harmonics. Spherical harmonic coefficients from partial sky coverage 
is related to the true ones as follows: a/ m = ^Ww mm ia}jf^,, where Ww mm i = 
J d 2 hW(h)Y l ^(ii)Yi' m i(h), and W(h) is zero in the masked region and one elsewhere. 
To estimate the amount of leakage, we assume CMB is purely Gaussian and has zero 
correlation: 

We find the expectation value of Df of the pure Gaussian CMB is 

(A 3 >-^tEK^) ( A 
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In the third equality, we have taken Cy out of the summation and equate it to C; + i, 
since W i p mm ,W / i + 2,/' m m' peaks sharply around V = I + 1, while Cp varies much slowly 
in comparison to Wl Vmm ,Wi+2,i'mm' ■ in the fourth equality, we have used the identity 
Ylii'm' ^'m'( n )^W'( n ') = <5(n — n'). From Eq. IA.14 we see the leakage from Gaussian 
CMB to Df can be estimated by: 
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The value as high as unity indicates that most of Gaussian CMB power is leaked into 
the fictitious correlation, while a zero value indicates no leakage. We have estimated 
Df leakages for the WMAP team's KQ75 and KQ85 mask and show them in Table 
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Table Al. D3 Leakage 
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IA11 We find Df leakage is as high as ~ 20% and other leakgages are ~ 2%, which 
are all greater than the expected correlations, given Alfven wave parameters in Table 
[TJ Therefore, we find CMB maps, which requires the KQ75 or KQ 85 mask, are not 
suitable for estimation of correlations induced by cosmological Alfven waves. 



